load(file = "Rdata/frap_2020_dis_Rx.Rdata")
load(file = "Rdata/CalMapper_activities_fire_2020.Rdata")
load(file = "Rdata/frap_2020_Rx.Rdata")
load(file = "Rdata/NF.Rdata")
It’s too impossible to actually find overlaps of activity acres so I’m going to stick with dissolving them and finding footprint acreage overlap
act_20_diss <- act_2020 %>%
st_union()
act_20_diss <- st_transform(act_20_diss, st_crs(frap_20))
overlap <- st_intersection(frap_diss_2020, act_20_diss)
overlap_ha <- as.numeric(st_area(overlap)/10000) %>%
round()
overlap_ha
## [1] 6219
frap_20_ha <- as.numeric(st_area(frap_diss_2020)/10000) %>%
round()
calmapper_20_ha <- as.numeric(st_area(act_20_diss)/10000) %>%
round()
paste("total area of FRAP dissolved polygons is", frap_20_ha, "ha")
## [1] "total area of FRAP dissolved polygons is 12100 ha"
paste("total area of CAL MAPPER dissolved polygons is", calmapper_20_ha, "ha")
## [1] "total area of CAL MAPPER dissolved polygons is 17478 ha"
mapview::mapview(list(act_20_diss, NF, frap_diss_2020),
col.regions=list("red","blue", "green"),
col=list("red","blue", "green"),
alpha.regions = 0.4)